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ABSTRACT 

New results on the anisotropy of the cosmic microwave background (CMB) and its polarization based upon 
the first three years of data from the Wilkinson Microwave Anisotropy Probe (WMAP) have revised the electron 
scattering optical depth downward from Tes = 0. 17;!;jj q^ to Tes = 0.09 ± 0.03. This implies a shift of the effective 
reionization redshift from z^ — 17 to Zr ~ 11. Previous attempts to explain the high redshift of reionization 
inferred from the WMAP 1-year data have led to widespread speculation that the sources of reionization must 
have been much more efficient than those associated with the star formation observed at low redshift. This 
is consistent, for example, with the suggestion that early star formation involved massive. Pop III stars which 
early-on produced most of the ionizing radiation escaping from halos. It is, therefore, tempting to interpret the 
new WMAP results as implying that we can now relax those previous high demands on the efficiency of the 
sources of reionization and perhaps even turn the argument around as evidence against such high efficiency. We 
show that this is not the case, however. The new WMAP results also find that the primordial density fluctuation 
power spectrum has a lower amplitude, cg, and departs substantially from the scale-invariant spectrum. We 
show that these effects combine to cancel the impact of the later reionization implied by the new value of Te, on 
the required ionizing efficiency per collapsed baryon. The delay of reionization is surprisingly well-matched 
by a comparable delay (by a factor of ^ 1 .4 in scale factor) in the formation of the halos responsible for 
reionization. 

Subject headings: cosmic microwave background - cosmology: theory - diffuse radiation - galaxies: formation 
- intergalactic medium 



1. INTRODUCTION 

One of the most important outstanding problems in cosmo- 
logical structure formation is how and when the universe was 
reionized. Observational constraints such as the Thomson 
scattering optical depth to the last scattering surface (Kogut 
et al. 2003; Page et al. 2006) from the large-angle polariza- 
tion anisotropy in the CMB detected by WMAP and the in- 
tergalactic, hydrogen Lya absorption spectra of high-redshift 
quasars (e.g., Becker et al. 2001) provide crucial constraints 
on the theory of cosmic reionization and the structure forma- 
tion which caused it during the early epochs that have thus 
far escaped direct observation. The WMAP first-year data im- 
plied an electron scattering optical depth, Tgs = 0.17, which 
seemed surprisingly large at the time, since it was well in 
excess of the value, Tes — 0.04, for an intergalactic medium 
(IGM) abrubtly ionized at Zr — 6.5, the reionization epoch 
which had been suggested by quasar measurements of the 
Gunn-Peterson (Gunn & Peterson 1965; "GP") effect. In or- 
der for such an abrupt reionization to explain the high value 
of 0.17 observed by WMAP for Tes, in fact, Zr — 17 is re- 
quired. This presented a puzzle for the theory of reionization: 
How was reionization so advanced, so early in our observed 
ACDM universe, and yet so extended in time as to accumu- 
late the high Tgs observed by WMAP, while ending as late as 
z ~ 6.5 to satisfy the quasar spectral constraints? 

This stimulated widespread speculation regarding the effi- 
ciency for the formation of the early stars and/or miniquasars 
which were the sources of reionization, as well as for the es- 
cape of their ionizing photons into the IGM (e.g., Haiman 
& Holder 2003; Cen 2003; Wyithe & Loeb 2003; KapHng- 
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hat et al. 2003; Sokasian et al. 2004; Ciardi, Ferrara, & 
White 2003; Ricotti & Ostriker 2004). A general consen- 
sus emerged that the efficiencies for photon production and 
escape associated with present-day star formation were not 
adequate to explain the early reionization implied by the high 
Tes value, given the rate of early structure formation expected 
in the ACDM universe. Common to most attempts to explain 
the high Tes was the assumption that early star formation fa- 
vored massive Population III stars, either in "minihalos," with 
virial temperatures Tvii < 10'* K, requiring that H2 molecules 
cool the gas to enable star formation (Abel, Bryan, & Nor- 
man 2002; Bromm, Coppi, & Larson 2002), or else in larger 
halos with T^i^ > 10"* K, for which atomic hydrogen cooling is 
possible, instead. A high efficiency for turning halo baryons 
into stars and a high escape fraction for the ionizing radiation 
into the IGM were generally required as well. Several effects 
were suggested that could extend the reionization epoch, too, 
including the rising impact of small-scale structure as a sink 
of ionizing photons (e.g., Shapiro, Iliev, & Raga 2004; Iliev, 
Shapiro, & Raga 2005; Iliev, Scannapieco, & Shapiro 2005), 
the suppression of low-mass source-halo formation inside the 
growing intergalactic H II regions (e.g., Haiman & Holder 
2003), and a general decline of the efficiency for releasing 
ionizing radiation over time (e.g., Cen 2003; Choudhury & 
Fen-ara 2005). 

With three years of polarization data, WMAP (henceforth, 
"WMAP3"), has now produced a more accurate determi- 
nation of Tes, which revises the optical depth downward to 
Tes = 0.09 ± 0.03 (Page et al. 2006). This value is consis- 
tent with an abrupt reionization at Zr = H, significantly later 
than that implied by the WMAP first-year data (henceforth, 
"WMAPl"'). It is natural to wonder if this implies that the 
high efficiency demanded of ionizing photon production by 
WMAPl, described above, can now be reduced, accordingly, 
to accommodate the later epoch of reionization determined by 
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Fig. 1 . — Halo abundance vs. mass for new (WMAP3) and old (WMAPl) 
parameters at z = 15, as labelled. Top: Press-Schechter mass function. Bot- 
tom: Fraction of matter /coii(> ^) in halos with mass greater than M. 

WMAP3. In what follows, we will show that this is not the 
case. 

Structure formation in the ACDM universe with the pri- 
mordial density fluctuation power spectrum measured by 
WMAP3 is delayed relative to that in the WMAPl universe, 
especially on the small-scales responsible for the sources of 
reionization. This, by itself, is not surprising, since there was 
always a degeneracy inherent in measuring the amplitude of 
the primordial density fluctuations using the CMB tempera- 
ture anisotropy alone, resulting from the unknown value of 
Tes. Higher values of Tes, that is, imply higher amplitude den- 
sity fluctuations to produce the same level of CMB anisotropy. 
This degeneracy is broken by the independent measurement 
of Tes made possible by detecting the polarization anisotropy, 
as well. Hence, when WMAP3 revised the value of Tgs down- 
ward relative to WMAPl, so it revised downward the ampli- 
tude of the density fluctuations. This same decrease of Tes 
implies a tilt away from the scale-invariant power spectrum, 
P(k) ex k"' with Hi = 1, which lowers the density fluctuation 
amplitude on small scales more than on large scales. As we 
will show, this delays the structure formation which controls 
reionization by just the right amount such that, if reioniza- 
tion efficiencies were large enough to make reionization early 
and Tes = 0.17 in the WMAPl universe, the same efficiencies 
will cause reionization to be later in the WMAP3 universe and 
Tes ~ 0.09, as required. 

In §2, we compare the rate of structure formation in ACDM 
according to WMAPl and WMAP3, on the scales relevant to 
reionization. In §3, we relate the history of reionization to the 
growth of the mass fraction collapsed into source halos, and 
use this to compare the reionization histories in WMAP3 and 
WMAPl universes. Our conclusions are summarized in §4. 

We adopt cosmological parameters {Vt,nh^, ilbli^, h, n,, 
CTg) = (0.14,0.024,0.72,0.99,0.9) (Spergel et al. 2003) and 
(0.127,0.022,0.73,0.95,0.74) (Spergel et al. 2006) for 
WMAPl and WMAP3, respectively. The most notable 
changes from old to new are: a reduction of normalization 
of the power spectrum on large scales (ct^ = 0.9 -^ 0.74) and 
more "tilt" (n, = 0.99 — > 0.95). Throughout this paper, we use 
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Fig. 2. — Collapsed fraction vs. redshift for halos with virial temperatures 
greater than r^in = 10* K, for WMAPl and WMAP3. as labelled. Dashed 
curve, nearly on top of "WMAP3" curve, is "WMAPl" curve with i+z ^ 

(l+z)/1.4. 

the transfer function of Eisenstein & Hu (1999). 

2. STRUCTURE FORMATION AT HIGH REDSHIFT 

A fundamental building block of models of reionization is 
the fraction of the mass contained in virialized halos - the 
"collapsed fraction" - the sites of ionizing photon production 
and release. Using the Press-Schechter formalism (Press & 
Schechter 1974), this collapsed fraction is given by: 



/coll (Z) = erfc I'min (z) I V2 



(1) 



where i'min(-z) = ^c/[Diz)(yiMj„i„)], a^{M) is the variance in 
the present-day matter density field according to linear pertur- 
bation theory, as filtered on the mass scale M, D{z) is the linear 
growth factor (Z)(z) oc 1/(1 +z) and Sc = 1.686 in the matter- 
dominated era), and Mmin(z) is the minimum mass for col- 
lapsed objects. For studies of reionization, the minimum mass 
is typically parameterized in terms of the minimum virial tem- 
perature, rnjin, of halos capable of hosting ionizing sources, 

M,™„~4xlO^M0[(7;„in/lO%)(lO/(l+z))(1.22/M)]^/^ 
where /i = 1 .22 for fully neutral gas (lliev & Shapiro 2001). 

For ACDM, a(M) for M - 10*^ - IO^Mq is lower for 
WMAP3 than for WMAPl by about 30 percent. During 
reionization, when such halos are still rare, we expect their 
abundance to be exponentially suppressed by this factor. This 
is clearly shown in Figure 1, where the new halo abundance 
and collapsed fraction are lower than the old ones by 1-2 
orders of magnitude. Since the threshold for halo collapse 
scales at these redshifts as 6c/D(z) oc 1 +z, structure forma- 
tion on these mass scales is delayed by a factor 1/0.7 ~ 1.4 
in scale factor. This is illustrated in Figure 2, where we plot 
/coll (7^ > Tmin = 10"^ K) vcrsus redshift and show that the shift 
of 1 .4 in scale factor provides an excellent description of the 
delay in structure formation which results. 

For the simplest possible reionization model, in which the 
universe is instantly and fully ionized at some redshift z, , the 
optical depth Tes oc (l + Zr)'^^^- If we assume that /coii(Zr) is a 
constant, so that reionization occurs when the collapsed frac- 
tion reaches some threshold value, then our simple estimate 
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implies that the change in the cosmological parameters alone 
reduces Tes by a factor of 1.4-'/^ = 1.65, from Tes ~ 0.17 to 
Tes ~ 0.1. In the next section we discuss the motivation be- 
hind tying the reionization history to the collapsed fraction. 

3. REIONIZATION fflSTORY 

An important quantity in the theory of cosmic reionization 
is the number of ionizing photons per hydrogen atom in the 
universe required to complete reionization^. In the absence 
of recombinations, this ratio is unity. Given some observa- 
tional constraint on the epoch of reionization, such as the on- 
set of the GP effect at z ~ 6.5, we can deduce that at least one 
ionizing photon per atom had to have been released by that 
time. This ratio can then be used to predict other quantities, 
such as the associated metal enrichment of the universe (e.g., 
Shapiro, Giroux, & Babul 1994) or the intensity of the near 
infrared background (e.g., Santos, Bromm, & Kamionkowski 
2002; Fernandez & Komatsu 2006). Most models of cosmic 
reionization link the ionized fraction of the IGM to the frac- 
tion of matter in collapsed objects capable of hosting stars 
(e.g., Shapiro, Giroux, & Babul 1994; Chiu & Ostriker 2000; 
Wyithe & Loeb 2003; Haiman & Holder 2003; Furlanetto, 
Zaldarriaga, & Hernquist 2004; Iliev et al. 2005; Alvarez et 
al. 2006). For a given model, reionization is complete when- 
ever the total number of ionizing photons emitted per hydro- 
gen atom reaches some threshold value. Along with the es- 
cape fraction, star formation efficiency, and stellar initial mass 
function, the evolution of the collapsed fraction /coii(z) forms 
the basis for calculation of this ratio and thus the reionization 
history. 

To relate Tes to the halo abundance encoded in /cou, it is 
necessary to determine the relationship between the reioniza- 
tion history and the collapsed fraction. If we assume every H 
atom which ends up in a collapsed halo releases on average 
f^(z) ionizing photons, and that e-^iz) is the number of ioniz- 
ing photons consumed per ionized H atom, then we can write 
a simple relation between /cou and the mean ionized fraction. 



Xe(z) = ^^/coll(z) = C(z)fco\\(z) 



(2) 



(e.g., Furlanetto, Zaldarriaga, & Hernquist 2004). For sim- 
plicity, we will assume a constant value, ((z) = Co (this sim- 
plification does not affect our main conclusions), and fix the 
value of Co for a given Tmin, so that Tes = 0.17 for WMAPl. 
Reionization is complete when the collapsed fraction reaches 
a threshold given by /coii(Z/-)Co = 1- 

In Figure 3, we plot the value of z/min which corresponds to 
Tmin = lO'* K. As mentioned in §2, WMAP3 implies a delay 
of structure formation by ^^ 1 .4 in scale factor. In the lower 
panel, we compare the reionization histories for WMAPl and 
WMAP3 according to equation (|2j, for the same efficiency Co- 
The same shift by a factor 1 .4 in scale factor is also present 
in the reionization histories, which is not surprising, since we 
have assumed that x^ ex /coii, and f^oii is a unique function of 
z/min- As mentioned in §2, this change can account for a shift 
in the implied value of Tes from 0. 17 to 0.1, quite close to the 
WMAP3 value of 0.09 ± 0.03. On the basis of this simple cal- 
culation, we conclude that the reduction of Tes from WMAPl 
to WMAP3 does not, itself, significantly reduce the demand 
for high efficiency of ionizing sources imposed previously by 
WMAPl. 

' For simplicity, we will neglect helium reionization. This does not effect 
our basic conclusions here. 
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Fig. 3.— Evolution with redshift for WMAPl and WMAP3, as labelled. 
Top: Threshold for collapse, I'min, for a halo with virial temperature 10"* K. 
Dashed curve, nearly on top of "WMAP3" curve, is "WMAPl" curve with 
1 +z ^ (1 +z)/1.4. Bottom: Reionization histories given by Xe = Co/coll. 
labelled by the corresponding values of Tcs, for fo = 170 and T^,^ = W* K 
(solid), and i^o = 35 and T^i^ = 2 X 10^ K (dotted). The two dashed curves are 
"WMAPl" curves with 1+z-^ (l+z)/1.4. 



3.1. Effect of recombinations 

Recombinations undoubtedly play an important role during 
reionization. To first approximation, they should determine by 
what amount the parameter e-y(z) appearing in equation Q ex- 
ceeds unity. The quantity e^- 1 is equal to the average number 
of recombinations that all ionized atoms must undergo during 
reionization, A',ec- As shown by Iliev et al. (2005), A',ec — 0.6 
at percolation for large scale simulations of reionization that 
resolve all sources with masses greater than w 2 x IO^Mq, but 
do not resolve clumping of the IGM on scales smaller than 
w 700 comoving kpc. Surely, smaller scale structure affects 
reionization strongly (e.g., Iliev, Shapiro, & Raga 2005; Iliev, 
Scannapieco, & Shapiro 2005), and therefore the number of 
recombinations per ionized atom is likely to be higher. For 
example, Alvarez, Bromm, & Shapiro (2006) found that the 
recombination time in the gas ionized by the end of the life- 
time of a 100 Mq star embedded in a IO^Mq halo at z = 20 is 
~ 20 Myr, roughly one tenth of the age of the universe at that 
time. 

At the high redshifts considered here, the ratio of the age 
of the universe to the recombination time is proportional to 
(1 +2)^/^. Since structure formation is later for WMAP3 than 
for WMAPl by a factor of 1.4 in scale factor, photon con- 
sumption due to recombinations is lower for WMAP3 by a 
factor ^ 1.4"'''^ = 1.65. Even if recombinations dominate the 
consumption of ionizing photons during reionization, there- 
fore, the new WMAP data require an efficiency Co which is at 
most a factor of only ^ 1 .65 lower than that for the first year 
data. This is true even if clumping increases toward lower 
redshift, since the evolution of clumping follows structure for- 
mation and is, therefore, similarly delayed. 
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4. DISCUSSION 

We have shown that the new cosmological parameters re- 
ported for WMAP3 imply that structure formation at high red- 
shift on the scale of the sources responsible for reionization 
was delayed relative to that implied by WMAPl. This de- 
lay can account for the new value in Tes without substantially 
changing the efficiency with which halos form stars. Recom- 
binations are fewer when reionization is later, but the reduc- 
tion is modest. Even the IGM clumping factor on which this 
recombination correction depends follows the delay in struc- 
ture formation. 

An important additional constraint on reionization is that 
it end at a redshift z > 6.5, in order to explain the lack of 
a GP trough in the spectra of quasars at z < 6.5. Because 
the GP trough saturates at a very small neutral fraction, the 
quasar data alone do not tell us when the universe became 
mostly ionized. Indeed, it is possible for the ionized frac- 
tion to have been quite high akeady at high redshift z ^ 10 
while there remained a neutral fraction sufficiently high to 
satisfy the GP constraint at z ~ 6.5 (e.g., Choudhury & Fer- 
rara 2006). While the universe may become mostly ionized 
well before z ^ 6.5, it cannot be later than this, however. Be- 
cause of the shifting in time of structure formation we have 
described, any model of reionization which previously satis- 
fied the WMAPl Tes ^ 0.17 constraint and became mostly 
ionized at z < 9 would now reionize too late to be compati- 
ble with the quasar observations. Recently, Haiman & Bryan 
(2006) used this fact to deduce that the formation of massive 
Pop III stars was suppressed in minihalos. 

While a^ = 0.744 and n, = 0.951 are the most likely values 
obtained from the new WMAP data alone, there remain sig- 



nificant uncertainties. When combined with other data sets 
such as large-scale structure (e.g., Spergel et al. 2006) and 
the Lyman-a forest (e.g., Viel, Haehnelt, & Lewis 2006; Sel- 
jak, Slosar, & McDonald 2006), important differences arise. 
While these differences may seem small from the point of 
view of the statistical error of the observational data, the im- 
plications for reionization can be quite dramatic, as we have 
seen here. It is also important to note that these measurements 
of the power spectrum are on scales much larger than those 
relevant to the sources of reionization. As such, the theory of 
reionization requires us to extrapolate the power spectrum by 
orders of magnitude beyond where it is currently measured. 
This means that the study of reionization is crucial to extend- 
ing the observational constraints on the origin of primordial 
density fluctuations (e.g., by inflation) over the widest range 
of wavenumbers accessible to measurement. Direct observa- 
tions of the high redshift universe such as 21 -cm tomography 
(e.g., Iliev et al. 2002; Zaldarriaga, Furlanetto, & Hernquist 
2004; Shapiro et al. 2006; Mellema et al. 2006) and large- 
aperture infrared telescopes such as JWST promise to dimin- 
ish the uncertainties which currently prevent us from making 
reliable statements about the nature of the first sources of ion- 
izing radiation. 
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